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ABSTRACT 

Physical properties of the star-forming regions in the local Luminous Compact 
Blue Galaxy (LCBG) NGC 7673 are studied in detail using 3D spectroscopic data 
taken with the PPAK integral field unit at the 3.5-m telescope in the Centro As- 
tronomico Hispano Aleman. We derive integrated and spatially resolved properties 
such as extinction, star formation rate (SFR) and metallicity for this galaxy. Our data 
show an extinction map with maximum values located at the position of the main 
clumps of star formation showing small spatial variations (E(B— V)t = 0.12— 0.21 mag). 
We derive a Ha-based SFR for this galaxy of 6.2 ± O.8M0/yr in agreement with the 
SFR derived from infrared and radio continuum fluxes. The star formation is located 
mainly in clumps A, B, C and F. Different properties measured in clump B makes 
this region peculiar. We find the highest Ra luminosity with a SFR surface densit y 
of 0.5 M0yr~^kpc~^ in this clump. In our previous work (|Perez-Gallego et al.ll2010[ ), 
the kinematic analysis for this galaxy shows an asymmetrical ionized gas velocity field 
with a kinematic decoupled component located at the position of clump B. This re- 
gion shows the absence of strong absorption features and the presence of a Wolf-Rayet 
stellar population indicating this is a young burst of massive stars. Furthermore, we 
estimate a gas metallicity of 12 + log(0/H) = 8.20 ± 0.15 (0.32 solar) for the inte- 
grated galaxy using the R23 index. The values derived for the different clumps with 
this method show small metallicity variations in this galaxy, with values in the range 
8.12 (for clump A) - 8.23 (for clump B) for 12 + log(0/H). The analysis of the emission 
line ratios discards the presence of any AGN activity or shocks as the ionization source 
in this galaxy. Between the possible mechanisms to explain the starburst activity in 
this galaxy, our 3D spectroscopic data support the scenario of an on-going interaction 
with the possibility for clump B to be the dwarf satellite galaxy. 

Key words: galaxies: starburst — galaxies: individual (NGC 7673) 



1 INTRODUCTION 

Starburst galaxies are those in which star formation 
and associa ted phenomena dominate the total energetics 
(IWeedmanl figsa). These galaxies have larger star formation 
rates (SFRs) per unit area than normal galaxies, and to 
produce their current stellar masses at their current SFRs 
they would take much less than the age of the Universe 
(|Kennicuttlll998bl ). Starburst galaxies are found at differ- 
ent redshifts, which denotes their cosmological relevance, 
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and turn those found nearby into perfect candidates to 
study the mysteries of the star formation process through- 
out time when they can be equally and properly selected 
at different epochs of the Universe . Both distant Lyman 
Brea k Galaxies fLBGs, [Steidel et al.lll996l : iLowenthal et al.l 
I997II and closer L uminous Compact Blue Galaxies (LCBGs, 
Werk et aLll2004l l fall into this category. 

LCBGs are, as described bv Werk et alj (|2004|V galaxies 
with (i) absolute blue magnitude (Mb) brighter than -18.5; 
(ii) effective surface brightness (SBe) brighter than 21 B- 
mag arcsec"^; and {iii){B — V) colour bluer than 0.6. In 
an observational parameter space defined by these observa- 
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tional properties (i.e., Mb, SBe, and {B — V)) LCBGs and 
distant LBGs share the same optical properties. Galaxies 
belonging to this population (i) are morphologically hetero- 
geneous, (ii) form stars at around 10 - 20 Mq yr~^, (iii) show 
velocity widths of 30 - 120 km s~^, (iv) are as compact as 
2-5 kpc, and (v) have metallicities lower than solar. 

LCBGs play an important role in galaxy evolution over 
cosmological time scales, as has been shown by various ob- 
servational studies in the last decade ( Werk et al. l20o3: 



Melbourne et al.1 l2007l : iGuzman etHI Il997l : ISteidel et all 
2QQ3I ). The lack of LCBGs in the local universe when com- 



pared to higher redshifts arises the issue of how this popu- 
lation has evolved in the last 9-10 Gyr. Two are the main 
suggested scenarios (not mutually exclusive): (i) LCBGs 
are the progenitors of today's spheroidal galaxies, low-mass 
(M < 10^° ^o) spheroidal galaxies or dwarf elliptical 
galaxies whose properties according to evolutionary mod- 
els would be matched by a typical LCBG after a 4 - 6 Gyr 
fading process showing low metalliciti es (|Koo et al.l 1 19941 : 
iGuzman et gUl 19981 : iNoeske et al.ll2006l ): and (ii) LCBGs are 
the progenitors of the spheroidal component of today's disk 
galaxies, present day small spirals, more massive (M ^ 10^° 
Mq) than inferred from virial masses , whose emission is 
mostly due to a vigorous central burst (|PhilliDS et al. I ll997l : 
iHammer et aDl200ll : IPuech et al.]l2006l ). This scenario pre- 
dicts higher metallicities than those likely to be observed in 
local dwarf galaxies ( Kobulnicky Zaritsk v 1999). 

Key ingredients on the discussion on whether LCBGs 
evolve one way or another is their mass and metallicity. A 
reliable determination of their masses and metallicities is 
necessary to properly place LCBGs within one evolution sce- 
nario or the other, or to understand what makes them evolve 
one way or another. Masses of LCBGs can be derived from 
their rotation curves and velocity widths, nevertheless one 
needs to be careful with those, since most of these objects' 
kinematics might not be coupled to their masses due to su- 
pernova galactic winds, and both minor and major mergers. 
The metallicity can be estimated in LCBGs by measuring 
the auroral emission line flux [OIII]A4363, which gives, along 
with the [OIII]AA4959, 5007 fluxes, a direct determination of 
the electron temperature. But beyond certain metallicities 
[OIII]A4363 detection is not possible and the oxygen abun- 
dance must be empirically determined using strong emission 
lines through the use, for example, of the R23 indicator. 

The knowledge not only of the integrated properties 
such as mass, metallicity, extinction or star- format ion rate, 
but their spatial distribution in the galaxy plays an impor- 
tant role to study the nature of LCBGs. Integral field spec- 
troscopy technique offers a great advantage to carry out this 
issue. In this way, we are able to derive full maps of physical 
properties such as extinction to locate the dusty regions in 
these galaxies. This is important to properly compute the 
extinction-corrected emission line maps and derive star for- 
mation rate and metallicity spatial distributions. 3D spec- 
troscopy makes it possible to carry out this work in an effi- 
cient way allowing us to test model predictions on the origin 
and evolution of these massive starbursts. 

Interest in LCBGs has multiplied following the initial 
observational results and the properties highlighted above, 
but their relation to today's galaxy population still remains 
unknown. In order to understand the nature of LCBGs and 
its role in galaxy evolution we have selected a representative 



sample of 22 LCBGs within 200 Mpc from the Sloan Digital 
Sky Survey (SDSS, lAdelman-McCarthv et aDl2006l\ Uni- 
versidad Complutense de Madrid fUCM. Zamor ano et aP 
1994) and Markarian catalogs (Markarian et al. 198^) that 
best resemble the properties of distant LCBGs, ensuring 
that this sample, although small, is representative of the 
LCBG population as a class by covering the whole range 
in luminosity, colo ur, su r face b rightness, and environment 
(see Perez- Gallego et al.1 (|201Q| ) (hereafter PGIO) for more 
details on the sample selection). We are carrying out a multi- 
wavelength study of this sample including not only the op- 
tical, but also the millimeter and centimeter ranges (e.g., 
ICarland et al.|[2007l V Arguably, one of the most important 
aspects of this study and the one we focus on this paper is 
the optical three dimensional (3D) spectroscopy. The optical 
is the best understood spectral range in nearby galaxies, and 
will be systematically studied at high- 2; with the new genera- 
tion of near-infrared multi-object spectrograph s and integral 
field units in 10-ni class telescopes (^Forster Schreiber et al.l 
I2OO6I : IPuech et al.1 120061 ). 

In this paper we focus on the physical properties of the 
star- forming regions of LCBG NGC 7673. We use integral 
field spectroscopy data which allow us to study both the 
integrated and the spatially resolved properties, such as ex- 
tinction, SFR and gas metallicity. The detailed analysis of 
these properties together with the kinematic analysis of the 
ionized gas discussed in PGIO, will help to understand the 
starburst nature of this galaxy. 



known as UCM2325+2318 and 
{z = 0. 011368) starburst galaxy 



NGC 7673, also 

MRK325, is a nearby ^ . 

widel y studied in the past f e.g.. [Dufiot-Augarde AlloinI 
19821: 
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Homeier Gallaghe [1999'). The small size, high surface 



brightness, strong emission lines and blue colours make NGC 
7673 a prototypical LCBG (see Table [1]). 

The starburst activity in this galaxy is located within 
a circle of 3.7 kpc radius (~ 2 x Re, where Re is the 
effective radius). It shows a clumpy structure with star- 
forming regio ns visible in optical as bri ght knots in the 
galactic disk (|Perez- Gonzalez et al ] I2OOI detected 87 dis- 
tinct star-forni i ng kn ots in a deep Ha image of this galaxy) . 
iHomeier et al.l (|200 2l identified 50 star cluster candidates 
working on HST/WFPC2 images, where the bluer and 
brighter clusters are strongly concentrated into the main 
starburst "clum ps" (see Figure [J) . From previo us ground- 
based images (^Dufiot-Augarde AlloinI Il982l ) SIX mam 
clumps have been identified A,B,C,D,E,F (see their Figure 
1) and will be referred to in this paper. 

The cosmology used in this paper is Ho = 70 km 
s"^ Mpc"\ The recession velocity of 3408 km/s from 
NASA/IPAC Extragalactic Database, implies a distance of 
d = 49 Mpc, and a projected scale of 250 pc per arcsec for 
this galaxy. This paper is structured as follows. Observations 
and data reduction are described in Section 2. Our measure- 
ments are shown in Section 3. The results and discussion are 
carried out in Section 4 and in Section 5. Finally, Section 6 
shows the summary of our work. 
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Table 1. NGC 767 3 Observational Properties. Redshift z from 
iHuchra et al Absolute B magnitude M(B), effective sur- 

face brightness SB e, colo ur index (B-V) and effective radius Re 
from lPisano e t al.' (EoOll) . 



Name NGC 7673 

z 0.011368 
M(B) -20.50 mag 

SBe 19.40 mag arcsec^ 
(B-V) 0.30 mag 

Re 1.9 kpc 



2 OBSERVATIONS AND DATA REDUCTION 

Objects from ou r sample were observed using PPAK 
(|Kelz et al.l l2006l l , a fibre based system for integral field 
spectroscopy operating at the 3.5-m CAHA telescope B The 
field of view of PPAK, as well as its good spatial and spec- 
tral resolution is optimum for the observation of our sample 
galaxies with typical half-light diameters of around 10 to 
30 , and velocity widths ranging from 60 to 120 km s~^. 
PPAK consists of 331 scientific fibres each of 2.7 in diame- 
ter, covering an hexagonal area of 74 x 65 on the sky. In 
addition, there are 15 calibration fibres and 36 fibres grouped 
in 6 bundles located at ^ 90 of the centre of the science 
bundle that are used to measure the sky simultaneously. 

PPAK observations of NGC 7673 were made over the 
nights of 2005 August 10 and 11 using two different setups. 
Firstly, a 300 lines mm~^ grating (V300) centred at 5316 
A was used. This low resolution configuration provided a 
spectral resolution of 10.7 A FWHM covering from 3600 to 
7000 A and allowed us to measure all the emission line ra- 
tios in one single spectrum, avoiding uncertainties associated 
to flux calibration and spectral response between different 
spectral ranges obtained with different configurations. Sec- 
ondly, a 1200 hues mm~^ grating (V1200) centred at 5040 
A was used. This intermediate resolution configuration pro- 
vided a nominal spectral resolution of 2.78 A FWHM (i.e. 
cr - 75 km s"^ at H/3), covering from 4900 to 5400 A. Three 
different dithering positions were observed in both spectral 
configurations (see Table 2 in PGIO for more observational 
details.) 

The data reduction of two-dimensional fibre spectra was 
done using the IRAf Q environment, R3D and Euro3D soft- 
ware ([Sanchez |2004| ) following the techniques described in 
ISanched (l2006h . The standard procedure consists on the 
following steps: bias subtraction, spectra extraction, wave- 
length calibration (see PGIO for more data reduction details 
in this step) , fibre to fibre correction and sky emission sub- 
traction. 

The aforementioned data reduction steps are done for 
each individual dithering exposure. The following steps in- 
stead are done using the information of the 3 dithering expo- 
sures to compute a final data cube absolute flux calibrated. 



^ Based on observations collected at the German-Spanish Astro- 
nomical Center, Calar Alto, jointly operated by the Max-Planck- 
Institut fiir Astronomic Heidelberg and the Instituto de As- 
trofisica de Andalucia (lAA/CSIC). 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories 



(i) Flux calibration. PPAK is based on pure fibre- 
bundles and does not cover the entire field-of-view, which 
imposes flux losses. Nevertheless our observational technique 
of dithering allows us to cover the complete field of view 
in three exposures. Therefore, our first approach is to de- 
termine a relative spectrophotometric calibration for each 
dithering exposure and recalibrate the spectra later, us- 
ing additional information coming from broad-band pho- 
tometry. The relative flux calibration requires the obser- 
vation of a standard spectrophotometric calibration star. 
We apply all the previous reduction steps to the calibra- 
tion star frame. We compare the measured spectrum of the 
star (flux coming from several fibres) with the absolute val- 
ues and determine the ratio between counts per second and 
flux. In this way we get the instrumental response which 
is applied to science frames to flux- calibrate them. To per- 
form an absolute flux calibration we compute broad-band 
PPAK maps in B, V and R filters. For B and R broad- 
band maps, fluxes obtained in different circular apertures 
over the PPAK field-of-view are compared with fluxes com- 
puted over JKT (B-band) a nd NOT ( R-band) Telescopes at 
La Palma, from Perez- Gonza lez et al . (2000, 2003) with the 
same apertures. JKT and NOT images were registered tak- 
ing as a reference the PPAK map and then convolved with a 
Gaussian filter t o match the spatial resolution of PPAK. For 
V-band filter the'Huchra (1977) tabulated value, my =13.11, 
is compared with the apparent magnitude computed from 
the PPAK V-band map. The comparison of relative PPAK 
fluxes with absolute fluxes in different photometric bands 
and in different apertures gives a mean factor of 0.44 with 
a standard deviation of 0.04, i.e PPAK fluxes need to be 
multiplied by this factor to match absolute fluxes. We com- 
pared the absolute PPAK Ha emission fluxes computed in 
different ap ertures with the fluxes from the Ha narrow-band 
image from lPerez-Gonzalez et al] (|2003l V In this case PPAK 
fluxes are ~ 15% smaller compared with those in the Ha 
narrow-band image. 

(ii) Locating the spectra in the sky. The location of 
the spectra in the sky is given by a position table that relates 
each spectrum with a certain fibre. Since we use the observa- 
tional technique of dithering exposures in order to cover the 
complete field-of-view, the next step in the data reduction 
is to create a data cube with the information for the three 
ditherings. This is done through the regular ization of the 
data: a grid of 1 arcsec/pixel is constructed over the field of 
view of the dithered exposures. Then, for each pixel of this 
grid, the spectra that covers the pixel are averaged (weighted 
by the fraction of pixel area in each dithering), obtaining a 
final data cube. This regular ization does not imply any in- 
terpolation. Finally, we obtain a data cube with a spectrum 
in each pixel where the exposure time for each pixel is differ- 
ent, depending on whether it has been covered by one, two 
or three fibres. 

(iii) Differential Atmospheric Refraction (DAR) 
correction. This correction is done by tracing the location 
of the intensity peak of a reference object in the field-of-view 
along the spectral range. These locations are estimated by 
determining the centroid of the object in the image slice 
extracted at each wavelength from the data cube. Then, 
the full data cube is shifted to a common reference by re- 
sampling and shifting each image slice at each wavelength 
using an interpolation scheme in the spatial direction. In the 
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10 arcsecs 





Fig;ure 1. Top: Ha im a^e from the NOT Telescope 
(IPerez-Gonzalez et al.l (|2003l )V Bottom: Optical image (WFPC2, 
F555W filter) from HST. The positions of the main clumps A to 
F are indicated. Both images are 40"x40" and the PPAK spaxel 
(2.7" diameter) is plotted at the right bottom corner. 

case of V300 configuration we apply maximum shifts of 0.2 
and 0.4 arcsec in the x and y axis, respectively. For V1200 
configuration we checked DAR correction is not required 
with shifts smaller than 0.1 arcsec. 



3 DATA MEASUREMENTS 

For each spectrum, different quantities such as emission line 
flux, continuum flux and equivalent width, were measured by 
fitting single Gaussian functions to the observed emission- 
line proflles using our own software. A two-component (emis- 
sion and absorption) gaussian flt was performed in order to 
correct the Balmer lines emission for underlying stellar ab- 
sorption (see Section [3?T]) . 

Here we analyze in detail the high resolution images 
obtained with the Hubble Space Telescope in the optical 
(WFPC2, F555W) and the NOT Telescope in Ha (see Fig- 
ure [1]) in order to match the spatial position of each clump 
with the starburst identifled in the Ha image. Clump A is 
composed by several distinct Ha knots, with three of them 
located in the central structure of the bar shape region. In 
the optical high resolution HST image two bright regions can 
be distinguished in this area. Clump B shows a strong Ha 



Figure 2. Top: PPAK B-band continuum map and Ho; emis- 
sion overlaid in contours. The positions of the main clumps 
A to F are indica t ed. B ottom: JKT B-band image from 
IPerez-Gonzalez et al.l (|2000l V North is up and East is left. 

emission located in one bright knot. The HST image shows 
clump B to be clumpy with a central peak of emission sur- 
rounded by two bright shells. It seems to be situated in a 
spiral arm connecting this region with the central structure. 
Clump C is another strong Ha emission region which shows 
several knots in the optical HST image. In Clump D the Ha 
emission is not seen at the cent re of the clump but at the 
edges. Clump E deflned in D uflot-Augarde AlloinI (|l982l l 
can be seen in the optical image as a region of faint emission 
connected to the central structure. Nevertheless this area is 
not associated with any Ha emission region in particular. At 
this location (SE of the galaxy) there are several knots of Ha 
emission. And flnally clump F is an extended star- format ion 
region com posed of many cluste rs and dominated by a sin- 
gle object (|Homeier et al.|[2002l ) embedded in a strong Ha 
emitting region. 

The spatial location of each clump in this paper is se- 
lected using the B-band and Ha PPAK maps (see Figure [2|). 
The spectrum corresponding to each clump is computed 
by adding up the spectrum from several pixels (see Fig- 
ure [3]). The integrated spectrum of the galaxy is computed 
co-adding the spectra of all pixels in PPAK dataset over the 
full fleld-of-view. 

The properties computed for different clumps in the fol- 
lowing sections are derived from the analysis of the spectrum 
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Figure 3. Upper panel: Spectra (in logarithmic scale and arbi- 
trarily scaled for better visualization) of different clumps in NGC 
7673 obtained summing up the spectra from some pixels defining 
these regions. Bottom panel: Integrated spectrum (in logarithmic 
scale) of NGC 7673 obtained co-adding all individual spectra of 
PPAK dataset over the full field-of-view. All the spectra shown 
here are obtained with the V300 grating configuration providing 
a nominal spectral resolution of 10.7 A FWHM. 



corresponding to each clump. Moreover, we have checked the 
agreement with the same properties derived from the mean 
values obtained from the computed map in the regions where 
the clumps are located. 



3.1 Stellar absorption correction 

To obtain an accurate value of the fluxes of the Balmer emis- 
sion lines we take into account the presence of an underlying 
stellar absorption. The stellar absorption can considerably 
affect the calculation of quantities such as the extinction. 
To solve this problem, two Gaussian functions are fitted in 
those cases where the absorption wings are visible (for ex- 
ample, that is the case for the H/3 Balmer line). For this 
purpose we use the V1200 configuration data cube whose 



spectral resolution (FWHM ^ 2. 5 A) allows us to better fit 
the absorption wings in the H/3 Balmer line. 

For each spectrum with enough signal-to-noise (SNR 
larger than 15 in the continuum), two different fits are per- 
formed. First, the H/3 Balmer line is fitted only with a Gaus- 
sian model in emission. Second, the same data is fitted with 
two Gaussian models, one in emission and another in absorp- 
tion. We tied the centroids of the gaussian in emission and 
absorption to be the same. The gaussian absorption width 
is upper limited to the value found when the total galaxy 
spectrum is fitted and pixels located at the wings of the H/3 
absorption are double weighted. The best fit is based in a 
scheme, and we accept the two gaussian fit as the best fit 
when the reduced improves at least a 15% with respect to 
the single gaussian fit. All the spectra in our data cube are 
fitted with two gaussian except the spectra corresponding 
to clump B, where we fit only a single gaussian in emission 
and then the H/3 absorption equivalent width is considered 
to be zero or negligible. 

The goodness of our fitting method is checked using 
simulated spectra. With signal-to- noise larger than 15 in 
the continuum, this method is able to measure absorption 
equivalent widths as small as ~ 2. 5 A with relative uncer- 
tainties smaller than 15%. For larger absorption equivalent 
widths the uncertainty decreases: for EWabs > 4 A we find a 
less than 10% relative error. 

In FigureUthe H/3 absorption equivalent width (EWabs) 
map and some representative fits in different regions are 
sho wn. The EWabs values are in the range - 6.6 A. 

iDufiot -Augarde Alloi 3 Cl982) employed a crude 
method to find EWabs ~ 8 A for the H/3 Balmer line in 
clump A and find no evidence of absorption wings in clumps 
B, C, and D. In Table [2] we summarize the H/3 absorption 
equivalent widths found for different clumps and the inte- 
grated galaxy spectrum. These values are measured fitting 
the spectra corresponding to each clump. While in clump 
B the underlying Balmer absorption is clearly less impor- 
tant, in other regions the stellar absorption is more signif- 
icant, with computed values of EWabs ~ 5.6 A in clump 
A, or 6.5 A A in regions between clumps A and D. When 
the integrated galaxy spectrum is fitted the H/3 absorption 
equivalent width is 4.8 di 0.5A. 

3.2 Emission Line Fluxes 

In this paper we focus on the study of the physi- 
cal properties of the gas through the analysis of the 
observed emission lines. In the V300 setup configura- 
tion the main emission lines observed are: [OII]A3727, 
H/3, [OIII]AA4959,5007, [HeI]A5876, [OI]A6300, Ha, 
[NII]AA6548,6584, [SII]AA6717,6731. The wide wavelength 
range of this configuration allows us to derive not only line 
emission fluxes maps but line emission ratios in a consistent 
way. 

The emission lines in each spectrum of the data cube 
were fitted by single gaussian profiles using our own soft- 
ware tools. We calculate the error in the line fluxes from 
the quadratic sum of the error coming from the absolute 
flux calibration and the error from the expression a\ — 
crcN^/2[l + EW/NA]^/2 (|Castellanodl200Ql ) where ac is the 
standard deviation in a box centred close to the measured 
emission line; N is the number of pixels used in the mea- 
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Figure 4. stellar absorption equivalent width map computed from the double gaussian fit method described in Section [3711 Contours 
from Hq; emission are over-plotted in white colour. Fitted spectra from different regions in NGC 7673 are shown. 



surement of the line flux; EW is the equivalent width of the 
line and A is the wavelength dispersion in A/pixel. Only the 
spectra with a minimum S/N of 15 in the flux detection are 
considered to ensure a relative uncertainty less than 10% in 
the flux determination. 

Some of the emission lines are fitted simultane- 
ously (Ha + [Nil] AA6548, 6584, [OIII]AA4959, 5007, 
[SII]AA6717,6731). In this case, the best fit is found 
by forcing all the lines fitted simultaneously to have the 
same velocity shift and width. The continuum under the 
emission line was fitted at the same time using a linear or 
quadratic function. 

We perform a pixel to pixel correction to the H/3 and 
Ha emission line fluxes computed from the V300 data cube, 
using the stellar absorption EW computed from H/3 Balmer 
line gaussian fits in the V1200 data cube as explained in 
section 13.11 We assume the stel lar absorption correction in 
H/3 and Ha are the same ([K^ucd (|l992l ) established that 
the Ha and H/3 equivalent widths are equal within a 30 
per cent uncertainty). In the case of Ha emission, the stel- 
lar absorption corrected fluxes are not significantly different 
from the uncorrected fluxes (relative differences are in the 
range 0-15% with 5% relative difference for the integrated 
galaxy flux). For H/3 emission fluxes the differences are much 
larger with relative differences over 30% for the integrated 
galaxy flux and the majority of pixels except for clump B 
(0%) and clump F (15%). To estimate the error in the Ha 
and H/3 absorption corrected fluxes we take into account the 
uncertainties from the absorption equivalent width and con- 
tinuum flux. The main source of error is coming from the 
uncertainty in the absolute flux calibration. 



3.3 Maps generation 



Finally we generate maps of spectral features such as emis- 
sion line fluxes, stellar absorption equivalent widths, ex- 
tinction and other properties from the regularized flux- 
calibrated data cube using our own software tools. Once we 
have a collection of values with the information of the mea- 
sured spectral feature in each pixel we can represent them 
using the relative position of each pixel in the sky. The com- 
puted maps consist of a 72 x 63 pixels grid where the spatial 
resolution element (fibre size of 2.7") is sampled by approx- 
imately 3x3 pixels. We decide to work with the spectral 
information in each pixel instead of in each fibre due to the 
application of DAR correction as it is explained in Section [2l 



To be able to compare our generated PPAK maps with 
other images we perform an astrometric calibration of our 
PPAK images using the Ha emission PPAK map and the Ha 
narrow-band image from NOT ( Perez- Gonzalez et al.ll2003h 
with known astrometry. The NOT image is convolved with 
a gaussian kernel function to match the PPAK spatial reso- 
lution. The positions of different HII regions in both images 
PPAK and NOT are used to compute an astrometric solu- 
tion using several tasks in IRAF package. The astrometric 
calibration is done by forcing the same x and y scale (1 arc- 
sec/pixel) with an error given by a rms of 0.09 arcsec and 
0.14 arcsec in the x and y axes respectively. 
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Table 2. Summary of PPAK results for NGC 7673 and the main clumps. Col. 1: Region from which the spectrum is analysed. Col. 2 : 
Measured stellar absorption equivalent width for H/3 Balmer line in A. Col. 3: Measured emission equivalent width for Ho; Balmer line 
with typical errors of O.SA. Col. 4: Colour excess obtained from the F(Hq;)/F(H^) ratio. Col. 5: Extinction corrected luminosity in 
units of lO^-'^erg/s. Col. 6: Extinction corrected Ha luminosity in units of lO^-'^erg/s . Col 7: Star formation rate in M0/yr computed 
from extinction corrected Ho; flux. Col. 8: Star formation rate surface density in M0/yr/kpc^. Col 9, 10 and 11: Emission line ratios 
computed from the extinction corrected fluxes. 



Region 

(1) 


EWabs 
(2) 


EW(Ha) 

(3) 


E(B-V) 

(4) 


L(H^) 

(5) 


L(Ha) 

(6) 


SFR 

(7) 


^SFR 
(8) 


[Nil] A6584 
Ha 

(9) 


[OIII] A5007 
(10) 


[SII]AA6716,6731 
Ha 

(11) 


NGC 7673 


4.8±0.5 


106 


0.17±0.03 


2.7±0.4 


7.8±0.9 


6.2±0.8 


0.08 


0.15 


1.6 


0.27 


Clump A 


5.6±0.6 


104 


0.19±0.03 


0.36±0.05 


l.OiO.l 


0.8±0.1 


0.49 


0.17 


1.6 


0.24 


Clump B 


0.0 


427 


0.21±0.01 


0.52±0.06 


1.5±0.1 


1.2±0.1 


0.52 


0.12 


1.9 


0.18 


Clump C 


3.1±0.5 


111 


0.16±0.03 


0.11±0.02 


0.33±0.04 


0.26±0.03 


0.23 


0.13 


1.7 


0.24 


Clump D 


5.1±0.5 


111 


0.12±0.03 


0.034±0.005 


O.lOiO.Ol 


0.08±0.01 


0.09 


0.14 


1.7 


0.27 


Clump E 


5.2±0.5 


113 


0.20±0.03 


0.032±0.005 


0.09±0.01 


0.07±0.01 


0.10 


0.17 


1.2 


0.31 


Clump F 


4.1±0.4 


171 


0.20±0.02 


0.17±0.02 


0.48±0.05 


0.38±0.04 


0.38 


0.14 


1.7 


0.24 



4 RESULTS 

4.1 Extinction map 

From the purpose of correcting for dust extinction the ob- 
served emission lines, we have derived the total colour ex- 
cess (E(B — V)t, i.e. the sum of the intrinsic and foreground 
Galactic reddening) from the Balmer line ratio Hq;/H/3. The 
departure from the measured ratio of the theoreti- 

cally expected value is produced by the extinction of light 
due to a screen of homogeneously distributed interstellar 
dust. 



E{B - V)t = Ko.p log 



2.86 



(1) 



where = [-0.4[k(Ha)-k(H/3)]]"^ is the ratio of 

the observed Ha and H/3 emission lines fluxes; and 2.86 rep- 
resents the expected Ha/H/3 value for case B recombination 
at T = lO^K and electron density of 100 cm~^ (jOsterbrockl 
fT989). Ka/3 depends on the selected extinction curve and 
takes the value 2.328 in CardelU et al. (1989). 

The Ha and H/3 stellar absorption corrected fluxes are 
used to derive the colour excess map (see Figure [5|). The fore- 
ground Galactic reddenin g at the position o f NGC 7673 is 
E(B - V)^^^ = 0.04 mag (Schlegel^eialj[l998i) • The PPAK 
colour excess map shows values in the range 0.04 - 0.35 mag 
in the majority of pixels. The dust extinction distribution 
seems to be related with the position of the main clumps of 
star formation (except for clump D) showing small spatial 
variations. In clumps A, B, E, and F the total colour excess 
is E(B - V)t ^ 0.20 ± 0.03 while for clumps C and D is 0.16 
and 0.12 respectively. In clump B we find relatively higher 
values for the colour excess (E(B — V)t ~ 0.4) in a region 
displaced with respect to the maximum Ha emission to the 
east direction. 

These valu es are in ag;reement w i th the colour excess 
determined by IPasauali Castangial (|2008l ) using cluster 
colours in different HST/WFPC2 bands, where the intrinsic 
cluster reddening E(B — V)i found is lower than 0.4 mag 
with the large majority of clusters in this galaxy having 
0. ^ E(B — V)i ^ 0.25 mag. They estimate cluster ages 
based on WFPC2 images and classified the clusters in three 
samples: young sample with ages < 2.5 Myr, intermediate 




Aa (arc sec) 



Figure 5. PPAK total colour excess map computed from the 
Balmer emission lines ratio F(Hq;)/F(H^) using Cardelli extinc- 
tion law. Hq; emission is shown in white contours. 



one with ages in the ragne 2.5-8 Myr and the old sample 
with ages > 8 Myr. Despite of the large age uncertainty, 
they point out that the stars clusters in NGC 7673 seem 
to get younger as their distance in the north-east direction 
from clump A increases. 

From the Balmer decrement we determine for the in- 
tegrated galaxy spectrum E(B — V)t = 0.17 mag. The ex- 
tinction Ana computed from this colour excess is 0.40 mag 
( Ay = 0.49) using the extinction curve from Carde ilTet al l 
(|1989|). 

The UV wavelength can provide us an independent at- 
tenuation estimation. We estimate the attenuation in the 
FUV and NUV bands using the ratio between the total dust 
emission in the range 1 — lOOO/xm (TIR) and the observed 
stellar emission in the FUV and NUV bands. The TIR lumi - 
nosity for this galaxy was derived by ICalzetti et aD (|2000l l 
from IRAS and ISO data. The FUV and N UV luminosi- 
ties were taken from I Gil de Paz et al.l (|2007l V By applying 
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Figure 6. Star formation rate density map in M0/yr /arcsec^ 
units, computed from reddening-corrected Ha flux PPAK map. 
SFR are computed in regions enclosed by white colour. 



the calibrations of ICortese et al. 1 (|2008l ) we find the atten- 
uation in the FUV and NUV to be Afuv = 1.08 mag and 
^Nuv = 0.86 mag, respectively. The slope of the UV spec- 
trum or, equivalent ly, the FUV— NUV colour, is known to 
be tightly correl ated with the UV extinction in starburs t 
galaxies (see e.g ICalzetti et al.1 Il994l : iMeurer et all 1 19991 ). 
Considering that NGC 7673 has FUV - NUV = 0.25, our 
values of Apuv and Anuv are fully consistent with what is 
expected for a starburst galaxy with that UV colour. 

We can use the previous values of the attenuation in 
the UV to estimate the attenuation in the V band. While in 
both cases we are referring to the extinction in the stellar 
continuum, stars contributing to the optical emission are 
usually less embedded within the dust than those domi- 
nating the UV emission. To take this fact into account, we 
combine the MW extinction curve of ICardelli et al.l (|l989l l 
wi th the simple sandwich model for dust geometry described 
in lBoselli et al.l ([2003), resulting in Ay = 0.35 mag. 

The extinction computed from the Balmer decrement 
is Av = 0.49 mag. This value corresponds to the extinction 
for the gas, it has t o be c orrected with the 0.44 factor found 
bv ICalzetti etHI (IT993 ) to be compared with the extinc- 
tion derived from the UV colours. After this correction, the 
stellar extinction results 0.13 mag lower compared with the 
attenuation in the V band using UV colours. Anyway, the 
0.44 factor is found for most dusty starbursts and this is not 
the case for NGC 7673 as it is shown from the low extinction 
values computed in this section. 



4.2 Star formation rates 

Different SFR measures based on Hq emission can b e found 
in the literature for this galaxy: ISchmitt et aP ([2006') derived 
a lower limit for the SFR of 1.3M0/yr for an observed Ha 
flux that was correcte d for Galactic extinct i on bu t not for 
intrinsic reddening. In^Pasau ali &: Castane^ial (|2008l l. the ob- 
served Ha flux from McQuad e et al.l ([l995|) is de-reddened 
using E(B-V)=0.54 mag (derived from the Balmer decre- 



ment using Ha and H/3 fluxes in iMcQuade et al ] ([19951 11 
and E(B-V)=0.36 mag (estimated for the exponentially 
decaying star formation history). They finally obtained a 
SFR between 3.3 and 5.2 Mrr. /vr. [Pisano et al.1 ([2001) de- 
rived a value of SFR=23.5M0/yr based on ICallego et al.l 
([ 19961 ) long-slit spectroscopic data. After reviewing how 
th is SFR value i s obt ained we find the L(Ha) published 
in Gallego et al. (1996) is highly overestimated due to the 
aperture c orrection applied. The same authors revisited this 
object in IPerez-Gonzalez et al.l ([2003l l and the SFR de- 
rived from the Ha imaging study is SFR=9 Mo/yr. Fur- 
ther more^_?TO^_7673_Jsalso^ well-known lum i nous FIR 
source (Sander s k Mirabell 1 1996V 'Ga rland et al.1 (|2005l ) de- 
rived a IR-SFR of 5.5 M© /yr (IRAS data) in agreement with 
the derived SFR from the 1.4 GHz r adio c ontinuum in NGC 
7673 using the prescription of^BdJ ([2003l V 

From the integrated galaxy spectrum we compute 
a total extinction-corrected luminosity L(Ha) = (7.8 it 
1.0) X 10^^ ergs" \ from which a SFR = (6.2 zb 0.8)Mc:)yr ~^ 
is derived (using the relation from Kennicutt (Il998al)) . 
The Ha luminosity derived bv Perez- Gonzalez et al.l ([2003l l 
(L(Ha) = 11.39 X 10^^ ergs" ^) is higher compared with our 
value mainly due to the absolute flux calibration (see Sec- 
tion [2|) and the extinction correction applied. They used the 
ratio F(Ha)/F(H/3) = 4.21 to correct the observed fluxes, 
and we find for the integrated galaxy spectrum the ratio 
F(Ha)/F(H/3) = 3.4. Our Ha-based SFR estimation is in 
agreement with t he SFR derived fro m infrared and radio 
continuum fluxes ([Garland et al.ll2005l . SFR = 5.5 Moyr"^). 

The SFR can also be e stimated fr o m the [OII]A3727 
luminosity. We use the iKennicuttI (19983) relation 
(SFR(M0yr-^) = (1.4=b0.4) x 10-^^L([OII])(ergs-^)), where 
the observed [OH] luminosity must be corrected for extinc- 
tion; in this case the extinction at Ha because of the manner 
in which the [OH] fluxes were calibrated to obtain that re- 
lation. We compute a total [OII]-based SFR of 8.5 Moyr"^, 
resulting in the range 6.1-11.0 Moyr"^ when we take into 
account the error associated in the relation L([OII])-SFR. 

We derive a SFR surface density map based on Ha 
extinction -corrected PPAK fluxes (see Figure O using the 
relation in IKennicuttI ([l998al l. The star formation is located 
mainly in clumps A, B, C and F as can be seen in Figure [6l 
Clump B shows the highest Ha luminosity with a SFR sur- 
face density of 0.52 M0yr~^kpc~^. Clumps D and E are re- 
gions with low Ha luminosity and their contribution to the 
total star formation rate is not significant. Although clump 
E was sel ected in this work for compariso n reasons with the 
results bv lDufiot-Augarde All^ ([l982l ) , this region is not 
associated with any star- forming region in particular, show- 
ing a low SFR. More than 50% of the total star formation 
rate in this galaxy is not located in the defined clumps. The 
star formation rates derived from the integrated galaxy spec- 
trum and the spectra corresponding to the different clumps 
marked in Figure [6] are tabulated in Table [21 



4.3 Emission line ratio maps 

We have computed different emission line ratios maps (i.e. 
[OIII]A5007/H/3, [NII]A6584/Ha, [SII]/Ha) to study the 
possible presence of an AGN or signs of shocks in this galaxy 
(see Figure O . The position of the different clumps except 
for clump E are very well traced in the [SII]/Ha map with 
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Figure 7. PPAK emission line ratio maps: a) [Nil] A6584/Hq;, b) 
[OIII]A5007/H^ and c) [SII]AA6716, 6731/Ha. The position of the 
main clumps are overlaid. 

clumps corresponding to regions of low emission line ratio 
(0.24). Clump B in this map shows the lowest value (0.18) 
and clump D the highest one (0.27). Region G is defined as 
a new region between clumps A and D. In this region the 
presence of the underlying stellar population is most notice- 
able (see Figure m of absorption stellar equivalent width) and 
the [SII]/Ha ratio reaches the value of 0.40. The map corre- 
sponding to the ratio [OIII] A5007/H/3 shows higher values, 
corresponding to higher ionization regions, at the position of 
clumps with values of 1.6, showing a slightly higher value 1.9 
at the position of clump B. The different clumps show values 
of [NII]A6584/Ha very similar in the range 0.12 (for clump 
B) and 0.17 (for clump A) . Only region G shows a higher 
value of 0.20 for this ratio. IPuflot-Augarde Allo"inl (|l982|) 
found a similar distribution of the [NII]AA6458, 6584/Ha ra- 
tio. In their analysis they computed the highest ratio value 
in clump A and lowest one in clump B. In Figure Elthe di- 
agnostic diagram [NII]A6584/Ha versus [OIII]A5007/H/3 is 
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Figure 8. Diagnostic diagram [Nil] A6584/Hq; versus 
[OIII 1A5007/H/3 is shown f or SDSS galaxies (grey pixels) 
from Brinchmann et al ] (120041) . The values for every pixel in the 
PPAK map are included in black crosses. The location of clumps 
A-F is also shown in an included figure. 
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Figure 9. Sum of the spectra in clump B where the WR stel- 
lar population is detected. The blue and red WR bumps can be 
observed. 



shown for every pixel in the PPAK map. For comparison, 
the values for SDSS galaxies are plotted in the same figure. 
We find values of [NII]A6584/Ha ratio in the range 0.1-0.2 
and [OIII]A5007/H/3 values in the range 1-2.5. The location 
of clumps A-F in this diagram is shown in the same Fig- 
ure and situate the clumps in a sequence of metallicity (as 
suggested by the se veral sets of theoretical photoionization 
models described in iDopita EvansI (|l986h l where clumps 
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Table 3. Metallicity estimations for each region derived from 
the R23 index. Col. 1: name of each region. Col. 2: R23 = 
log(([OII]A3727+ [OIII]AA4959, 5007)/H/3) values with typical er- 
ror of 0.07 dex. Col . 3: Derived oxyg en metallicity using R23 
method described in lKewlev Dopital (|2002l ) for the low metal- 
licity region, taking into account the ionization parameter. Col. 4: 
Measured Mb (corrected for Galactic extinction) in magnitudes 
with typical error of 0.1 mag. 



Region 


R23 


12+log(0/H) 


Mb 


NGC 7673 


0.72 


8.20±0.15 


-20.2 


Clump A 


0.68 


8.12±0.12 


-17.9 


Clump B 


0.73 


8.23±0.15 


-17.3 


B Surroundings 


0.76 


8.30±0.20 




Clump C 


0.73 


8.21±0.16 


-16.9 


Clump D 


0.69 


8.14±0.12 


-15.6 


Clump E 


0.69 


8.15±0.14 


-15.3 


Clump F 


0.73 


8.21±0.16 


-16.7 


Region G 


0.64 


8.07±0.11 





£0 and A would have higher metallicity values, clumps C, 
D and F intermediate metallicity values, and clump B the 
lowest one. These ratio values (see Table [2]) are clearly in 
agreement with those of HII regions of i ntermediate meta l- 
licity and with no signs of AGN activity (IOsterbrocklll989l ). 
The [SII]/Ha map is also computed showing values < 0.4 for 
all the pixels. The location of these values in the diagnostic 
diagram [SII]/Ha versus [NII]A6584 /Ha is compatible with 
a starburst galaxy (see Figure 12 in 'Rickes et all 1200 8') rul- 
ing out the presence of LINER or shocks as the ionization 
source in this galaxy. 

4.4 Wolf-Rayet stellar population 

Wolf- Ray et (WR) stars are evolved, massive stars, which 
are lo sing mass rapidly by means of a very strong stellar 
wind ('Maeder & Conti' 1994). The presence of WR stars can 
be inferred from the detection of WR bumps around A4650 
A (blue bump) and A5808 A (red bump), which are generally 
a blend of Hell and several metal lines. 

We detect the presence of WR stellar population in 
NGC 7673. The blue and red WR bumps can be observed 
in several fibres at the position of clump B. The sum of the 
spectra where WR bumps are detected is shown in Figure [9l 
The WR bumps intensity maximum coincides spatially with 
the position of Ha peak emission in clump B. The detection 
of a WR stellar population at this clump indicates this is a 
young (less than 4 Myr) burst of massive stars. I n agreement 
with this detection IPasquali Castangial (|2008h found clus- 
ter ages younger than 2.5 Myr and in the range 2.5-8 Myr 
for the majority of star clusters located in clump B. 

4.5 Metallicity 

iDuflot -Augarde AlloinI (|l982l l determined almost a con- 

^ Actually, clump E, although defined in 

iDufiot- Augarde AlloinI (|l982l V is not associated with any 
Ha emission region in particular. 



Table 4. Metallicity estimations for each fibre in Figure^] de- 
rived from the Te-method and R23 ratio as explained in Sec- 
tion ^5] Col. 1: fibre number. Col. 2: Oxygen metallicity com- 
puted from Te-method. Col. 3: Derived oxygen metallicity using 
R23 method. 



Fibre 


12 + l0g(0/H)Te 


12 + l0g(0/H)R23 


1 


7.96±0.05 


8.14±0.14 


2 


8.00±0.05 


8.14±0.14 


3 


> 8.02 


8.33±0.25 


4 


8.07±0.04 


8.06±0.11 


5 


7.95±0.06 


8.35±0.27 


6 


> 7.82 


8.31±0.22 


7 


> 8.01 


8.11±0.13 


8 


> 7.95 


8.28±0.19 


9 


8.10±0.05 


8.07±0.12 


10 


> 7.84 


8.12±0.13 



stant oxygen abundance in the clumps studied in 
NGC 7673 (12+lo g(0/H)^8.6) usin^ empi rical meth- 
ods d e scribed by Alloin et al.l (1 19791 ) and IPagel et"al] 
(|l979l ). IStorchi -Bergmann et al.l (Il994f ) obtained a simi- 
lar value of oxygen abundanc e (12+log(0 /H)=8.48) using 
the calibration diagrams from IPagel et al.l (jl979|) to obtain 
T[OIII]. 

A precise measurement of the weak auroral forbidden 
emission line [OIII]A4363 is needed to give an accurate de- 
termination of oxygen abundance in gaseous ionized nebu- 
lae. This emission line which is tem perature sensitive cor- 
relates with the overall abundance (IOsterbrocklll989l ). be- 
ing relatively strong in very low metallicity systems but 
undetectable for even moderately low metallicity galaxies 
(12 + log(0/H) > 8.3) (Denicolo et al. 2002). 

We inspected all the spectra in NGC 7673 and for 
some fibres the detection of the emission line [OIII]A4363 
is possible with enough signal-to- noise ratio (5cr-level). In 
particular we are able to detect [OIII]A4363 in a few spec- 
tra corresponding to clump B where the stellar absorption 
over the H7 emission is not so important and the emission 
flux for [OIII]A4363 can be measured reliably. We apply 
the Te-method to measure gas metallicity, using the rou- 
tine abund of the nebular package of the STSDAS/IRAF 
for this task. We assume a two-zone model for the star- 
forming nebula, a medium-temperature zone where oxygen 
is doubly ionized and a low-temperature zone where oxygen 
is assumed to be singly ionized and neutral. The electron 
temperature in the medium-temperature zone, Te(OIII) is 
derived fr om the [OII I](A4959+A5QQ7)/A4363 line inten- 
sity ratio (jlzotov et al.l ll994V The temperature in the low- 
temperature zone, Te(OII) is estimated using the relation 
te(OII) = 2[te(0III)"^ + 0.8]-\ proposed by Pagel eTell 
((l99i) based on HII region models of Stasinska (1990), where 
te are temperatures measured in units of lO'^K. The elec- 
tron density (Ne) is estimated from the [SII]A6716/A6731 
ratio. The computed total oxygen abundance is the sum of 
all measured oxygen ionic abundance, considering negligible 
the contribution of O^^ to the total abundance. 

The error estimation for the oxygen abundance comes 
mainly from the uncertainty in the [OIII]A4363 emission 
flux. We estimate 10-18% relative error in this emission flux. 
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Figure 10. Left: R23 = log(([OII] A3727 + [OIII] AA4959, 5007)/H/3) map with regions where the oxygen metaUicity (using R23 index) is 
computed overlaid. Right: Detail of the Ha emission in clumps B and C. The position of some fibres with the computed oxygen metallicity 
is overlaid. Oxygen metallicity (12+log(0/H)) is derived in five fibres located at clump B and C where the emission line [OIII]A4363 is 
detected (5-sigma level). In some fibres we compute a metallicity lower limit. 



We compute new electronic temperatures taking into ac- 
count this uncertainty, and derive new oxygen abundances. 
Relative errors in temperature range from 4% to 7% and 
result in final metallicity errors of 0.04-0.07 dex. 

In some fibres located at clump B, we use the signal-to- 
noise ratio of each spectrum to find an upper limit to the 
[OIII]A4363 line strength in the case this auroral line is not 
detected. With this estimation we can derive an upper limit 
to the electron temperature and hence a lower limit for the 
oxygen abundance. 

The oxygen abundance computed in clump B is in 
the range 12 + log(0/H)=7.96 - 8.10, with a value of 
12 + log(0/H) = 8.07 zb 0.04 solaiEl) for the fibre at 

the position of the maximum Ha intensity. Furthermore, the 
oxygen abundance derived from the sum of spectra in clump 
B is 12 + log(0/H) = 8.02 ± 0.05. We attempted to find an 
estimation of the metallicity summing up the signal from 
different fibres around clump B. In this case the [OIII]A4363 
line strength is difficult to measure due to the position of this 
emission line near to the H7 absorption stellar wing. Never- 
theless we estimate a lower limit for the oxygen abundance 
in the surroundings of clump B of 12 + log(0/H) > 7.87. 

A metallicity estimation for the rest of regions in 
this galaxy is obtained using the R23 index (R23 = 
log(([OII]A3727 + [OIII] AA4 959, 5007) /H/3). W e apply the 
R23 method described in iKewlev k Dopital (|20oj ) for 
the low metallicity region, taking into account the ion- 
ization parameter. From the spectrum corresponding to 
each clump we derived an initial guess for the metallic- 
ity u sing the [OH] /[OIII] and [NII]/[SII] emission Une ra- 
tios (|Charlot k, Longhettil l200lh . Taking a metallicity esti- 
mation for each clump, the first estimation of the ioniza- 

4 We consider [12+logrO/H)]rT.=8.69 lAllende Prieto et al.l(l200ll ) 



tion parameter is derived from the [OIII] /[OH] ratio. The 
R23 index is then used to determine the oxygen abundance 
and if this value is significantly different from the initial 
metallicity guess, the ionization parameter is again com- 
puted with the new metallicity estimation. This process is 
repeated until there is no significant change in the metallic- 
ity value which usually happens after the first iteration. The 
R23 values and the oxygen abundance estimations for dif- 
ferent regions in NGC 7673 are shown in Table [3] Figure [TOl 
shows the R23 map with values in the range 0.6-0.8. For the 
integrated galaxy spectrum we derive R23 = 0.72 ±0.07 and 
q = 2 X lO^cms-^ which giv es 12 + log(0/H) = 8.20 ± 0.15 
lower than previous works. Dufiot-Augarde AlloinI (|l982l l 
found oxygen abundance of 12 + log(0/H) = 8.45 for clump 
A and 12 + log(0/H) = 8.60 for clump B using the empir- 
ical method based on the ([OH] + [0III])/H/3 line intensity 
ratio and the opposite values when the [OIII] /[Nil] ratio is 
employed. They considered the two empirical methods pro- 
vide abundances in good agreement taking into account the 
uncertainty inherent to this kind of analysis. 

For comparison we computed the R23 oxygen abun- 
dance in those fibres where we detect [OIII]A4363 (see Fig- 
ure[TO]and TableH]). The R23 method overestimates the oxy- 
gen abundance in fibres 1,2, and 5 (differences of 0.14-0.40 
dex). The two methods give similar results in fibres 4 and 9 
where the signal-to- noise ratio in the [OIII]A4363 detection 
is higher suggesting that the R23 method is reliable. For the 
remaining fibres (3, 6, 7, 8 and 10) where the [OIII]A4363 
auroral line is not detected we only estimate lower limits for 
the abundances derived with the Te method. 

In conclusion, for clump A which dominates the galaxy 
morphology we derive a low R23 metallicity estimation (0.27 
solar). Clump B seems to be surrounded by a region which 
shows marginally higher metallicity (0.41 solar) compared 
with the metallicity derived in clump B (0.35 solar) from the 
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same method R23. Clumps D and E with the same metal- 
hcity estimation (0.28 solar) seem to belong to the same 
spiral arm. And finally clumps C and F situated at equidis- 
tant positions respect to clump A, have the same derived 
metallicities (0.33 solar). 



5 DISCUSSION 



5.1 Nucleus 



The loc ation of the nucleus in NG C 7673 is not well deter- 
mined. Homeier &; Gallagherl |l999) refers to clump A as the 
"nucleus" . If this is true we would expect a specific behaviour 
in the kinematic, velocity width maps and diagnostic dia- 
grams. From the analysis of the Ka kinematic map PGIO 
find the position of the kinematic centre in between clumps 
C and A. This location corresponds to the position of a 
maximum found in the gas velocity width map and roughly 
coincides with the photometric and geometric centres of the 
galaxy. 

This latter finding does not, however, rule out the possi- 
bility for clump A to host the nucleus. Indeed, our computed 
extinction map shows at the position of clump A a concen- 
tration of extinction although with not very high values. 
Probably the nucleus of this galaxy is located in this clump 
but the optical light does not penetrate enough to reveal the 
optical counterpart of the nucleus. 

Radio continuum emission in t his gjalaxy mig:ht help 
to locate the position of the nucleus. ICondon &: YinI (|l99d ) 
shows VLA data (A- array at 1.46 GHz) mapped at 3" and 
1.5" angular resolution. These radio continuum images show 
an irregular, diffuse morphology characteristic of a starburst 
with the two brightest radio peaks matching up he brightest 
optical clumps A and B. But there is no radio peak emis- 
sion at the position of the kinematic centre found in PGIO. 
Furthermore in their high-resolution radio map there is no 
evidence of any compact source stronger than S=0.5 mJy or 
of radio jets powered by a nuclear AGN. 

When we represent clump A in several diagnostics dia- 
grams we find no sign of any AGN activity (see Figure [8]). 
Mid-IR spectra can shed further light on the AGN versus 
starburst nature of the galaxy. In particular, the emission 
features due to Polycyclic Aromatic Hydrocarbons (PAHs) 
are known to be much weaker in A GNs than in starburst 
galaxies (see, e.g., IWeedman Houckll2009i . and references 
therein). In the particular case of NGC 7673, mid-IR spec - 
tra from the Infrared Sp ectrofflaph (I RS jHouck et al.|[20o3 ) 
onboard Spitzer (| Werner et al.l [20o3 ) are publicy available 
in the Spitzer archive (pro gram ID 73, PI: J . R. Houck). 
We used the tool CUBISM (jSmith et aLll2007l l to build and 
inspect the corresponding data-cube. PAHs features are vis- 
ible across the whole galaxy, and are particularly bright in 
clump A, whose spectrum resembles the starburst templates 
presented in Weedman Houckl (2009). Moreover, we do 
not observe strong high-ionization lines of [Ne V] or [O IV] 
which, if present, would point towards an AGN as the pri- 
mary ionizing source. Therefore, the mid-IR data supports 
the starburst nature of this galaxy. 



5.2 The peculiarity of clump B 

Different properties measured in clump B makes this region 
peculiar. This clump with bluer colour (Mb = —17.3) and 
medium-high star formation rate density (O.SMo/yr/kpc^) 
represents a compact decoupled kinematic component stud- 
ied in detail by PGIO. The HST image of this galaxy shows 
B region to be clumpy with a central peak of emission sur- 
rounded by two bright shells. Its bubble morphology is prob- 
ably the result of the stellar winds associated with massive 
stars. Indeed, we detect the presence of Wolf-Rayet stellar 
population in several fibres in clump B pointing out the pres- 
ence of massive stars, where young (less than 2.5 Myr old) 
and in termediate age cl usters (between 2.5 and 8 Myr) are 
found (|Pasquali &; Cast angia 2008). 

While at the centre of clump B we find extinction val- 
ues similar to clumps A and F, relatively higher values for 
the colour excess (E(B — V)t ~ 0.4) are found in a re- 
gion displaced respect to the maximum Ha emission to the 
east direction. Our measurements of the equivalent width of 
the underlying stellar population absorption features (e.g., 
H/3, H7) show a peculiarity at the location of clump B. 
The strength of these features is noticeable throughout the 
galaxy but it is almost absent at the location of clump B 
(see Figure |4|. The strength of the equivalent width of these 
features account for the age of the underlying population. 
Strong lines are mainly due to a population of A class stars 
and correspond to systems a few hundred thousand years 
old. Weaker lines are typical of either younger or older pop- 
ulations. Nonetheless, clump B shows signs of differentiation, 
at least marginally, with respect to its environment from the 
point of view of both, the gas and the stars. 

From our metallicity study clump B is a region with 
low metallicity 1/4 solar) surrounded by a region with 
marginally higher metallicity (see Figure fTOj) . Furthermore, a 
detailed study using the R23 index, shows some fibres at the 
east side of the central fibre in clump B with higher metal- 
licity (~ 1.2 solar). This region is located at the position 
where PGIO find maximum values for the velocity width. 

The location of NGC 7673 in the luminosity-metallicity 
relation is in agreement with the trend o bserved for 
intermediate-z star-forming galaxies from 
("2009) and from intermediate-z LCBGs from 



Salzer et al 



Hovos et al 



(2005) (see PGIO). But the properties of clump B, being 
less luminous and with approximately the same metallic- 
ity, po sitioned this region clo ser to the dwarf irregulars 
trend ([Richer McCalll Il995l ) as already seen in PGIO. 
But this region could also be a nearby giant HII region 
following the H/3 lu minosity- velocity width correlation by 
iMelnick et all (|l987h . Nevertheless, as explained in PGIO, 
clump B is brighter and more massive than any of the nearby 
giant HII regions within their sample, suggesting that an in- 
falling dwarf galaxy rather than a giant HII region, might 
be the reason for the peculiarities found in this clump. 

In summary, our 3D data show clump B is a compact 
region kinematically decoupled with high extinction and star 
formation rate density and with luminosity and metallicity 
properties that suggest a dwarf galaxy or a giant HII region 
nature for this clump. 
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5.3 On the nature of NGC 7673 

The star- format ion activity in this LCBG galaxy is located 
in several clumps as it is shown in the star formation rate 
surface density map. In general, the origin of the star- 
burst activity could be explained by different mechanisms 
which include, i nteractions and mergers_ (ISchweizerl 1 19871 : 
I Jog Dad I1992I ). bar instabilities (|Shlosman et alJ Il990ll 
and kinematic fr om SNe and stellar winds ([Heckman et alJ 
ll99Ql ). Although iHomeier" GaUaghed excluded an 

active interaction of NGC 7673 with its companion NGC 
7677 as the starburst trigger, they cannot rule out a past 
interaction. Another possibility suggested by these authors 
is the capture of a dwarf companion in a minor merger as 
the main starburst triggering which would account for the 
presence of the different bursts and the characteristic mor- 
phology in this galaxy. 

The HST image of this galax}0 shows a clearly distorted 
morphology with an inner disturbed spiral structure. Clump 
A dominates the galaxy morphology and probably the nu- 
cleus of NGC 7673 is hidden by extinction at this position. 
The ring-like shape with the absence of Ha emission at the 
centre of clump D (as seen in the high spatial resolution 
map) shows this is an evolved star-forming region decou- 
pled from the spiral structure where region E is located. 
Our kinematic study in PC 10 shows an asymmetrical ion- 
ized gas velocity map, where a decoupled kinematic compo- 
nent is found at the position of clump B. This region results 
peculiar not only from the kinematic study but also from 
different physical properties derived in this paper. This re- 
gion, with the highest Ha emission equivalent width and 
star formation surface density, is composed of young (we 
detect the presence of WR stellar population) and interme- 
diate age clusters where the underlying stellar population 
shows the absence of strong absorption features. We find no 
evidence for neither AGN activity (see Section [4. 3 p nor SNe 
galactic winds (PGIO) in this kinematically decoupled com- 
ponent, and from the metallicity and luminosity derived in 
this work, clump B is in agreement with being an extremely 
giant HH region or an in-falling dwarf g a laxy (see PGIO). 

As stated in iHomeier G allaghed (|l999l l. NGC 7673 
and NGC 3310, present morphological and spectroscopic 
similarities in H H and H I data which would lead us to 
think in the same triggering scenario for both starburst 
galaxies. NGC 3310 has been largely classified as a system 
which is undergoing a minor merg er with a dw arf compan- 
ion (see for example Mulder & v an Driell 1 19961 ) . A similar 
scenario would account for the inner bursts in NGC 7673 
and the different morphological features, such as arcs, in 
its outer parts. The spectroscopic and kinematic proper- 
ties derived for clump B in this work and in PGIO would 
lead us to think that this decoupled kinematic component 
is the dwarf satellite galaxy in the minor merger hypothesis 
althou js^h this needs to be c o nfirm ed with further studies. 
While IHomeier Gallagher! (|l99^ ) talk about a past mi- 
nor merger to allow the outer disk to mostly recover, the 
similarities with NGC 3310, the position of the decoupled 
kinematic component in the galaxy disk and the young and 
intermediate age of its starburst seem to support an on-going 
interaction. 

^ http:/ /www. spacetelescope.org/images/html/heic0205a.html 



6 SUMMARY 

Physical properties of the star-forming regions in the local 
Luminous Compact Blue Galaxy (LCBG) NGC 7673 are 
studied in detail using three dimensional data taken with the 
PPAK integral field unit at the 3.5-m telescope in the Centro 
Astronomico Hispano Aleman. From previous ground-based 
images six main clumps have been identified in this galaxy 
being four of the clumps strong Ha sources (A,B,C and F). 
We derive integrated and spatially resolved properties such 
as extinction, star formation rate and metallicity for this 
galaxy and our results are the following. 

1. Our data shows an extinction map with maxi- 
mum values located at the position of the main clumps 
of star formation showing small spatial variations. The 
colour excess values for the different clumps are in the 
range E(B - V)t = 0.12 - 0.21mag (0.17 mag for the in- 
tegr ated galaxy spectrum) , in agreement with the results 
from IPasquali Castangial (2008). 

2. We derive a Ha-based SFR for this galaxy of 6.2 zb 
O.8M0/yr in agreement with the SFR derived fro m infrared 
and radio continuum fluxes (jCarland et al.l [20051 ). The star 
formation is located mainly in clumps A, B, C and F. Differ- 
ent properties measured in clump B makes it peculiar. We 
find for this burst the highest Ha luminosity with a SFR sur- 
face density of 0.5 M0yr~^kpc~^. This region shows the ab- 
sence of strong absorption features and the presence of Wolf- 
Rayet stellar population indicating this is a young burst of 
massive stars. The analysis of the emission line ratios discard 
the presence of any AGN activity or shocks as the ionization 
source in this clump or any other in this galaxy. 

3. Furthermore, we estimate a gas metallicity of 12 + 
log(0/H) = 8.20±0.15 (0.32 solar), for the integrated galaxy 
using the R23 index. The values derived for the different 
clumps with this method show a small metallicity variations 
in this galaxy, with values in the range 8.12 (for clump A) - 
8.23 (for clump B) for 12 + log(0/H). 

The detection (5cr-level) of the auroral forbidden emis- 
sion line [OHI]A4363 in some fibres located in clump B al- 
lows us to determine metallicities using the Te-method and 
compare these values with the computed R23 oxygen abun- 
dance for the same fibres. Although in some cases the R23 
index overestimates the oxygen abundance (differences of 
0.14-0.40 dex), in those fibres with higher signal-to- noise ra- 
tio in the [OHI]A4363 detection both methods give similar 
results suggesting the reliability of the R23 method. 

4. Between the possible mechanisms to explain the star- 
burst activity in this galaxy, our 3D spectroscopic data sup- 
port the scenario in which a minor merger is taking place 
with the possibility for clump B, a compact decoupled kine- 
matic component, to be a dwarf satellite galaxy. 
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